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Abstract. We consider the 3-D Navier-Stokes initial value problem, 

Vt - uAv = -Vlv- \7v] + / , v{x, 0) = vo{x), X (*) 
where V is the Hodge projection to divergence-free vector fields. We assume 
that the Fourier transform norms ||/||;i(23) ^^'^ ll^o|lii(z3) ^.re finite. Using 
an inverse Laplace transform approach, similar to the earlier work we 
prove that an integral equation equivalent to (*) has a unique solution U{k, q), 
exponentially bounded for 5 in a sector centered on where q is the inverse 
Laplace dual to l/t" for n > 1. 

This implies in particular local existence of a classical solution to (*) for 
t G (0,T), where T depends on ||iio|lii and Global existence of the 

solution to NS follows if q) Hji has subexponential bounds as g — > 00. 

If / = 0, then the converse is also true: if NS has global solution, then 
there exists n > 1 for which necessarily decays. More generally, if 

the exponential growth rate in g of 1/ is a, then a classical solution to NS 
exists for t e (O, Q-^/"). 

We show that a can be better estimated based on the values of (7 on a 
finite interval [0, go]- We also show how the integral equation can be solved 
numerically with controlled errors. 

Preliminary numerical calculations of the integral equation over a mod- 
est [0, 10] q-interval for n = 2 corresponding to Kida ( 1201 ) initial conditions, 
though far from being optimized or rigorously controlled, suggest that this ap- 
proach gives an existence time for 3-D Navier-Stokes that substantially exceeds 
classical estimate. 

1. Introduction 
We consider the 3-D Navier-Stokes (NS) initial value problem 

(1.1) vt-iyAv = -V[v-\7v] + f{x), v{x,0) ^ voix), a; G ¥^[0, 27r], teR+ 

where v is the fluid velocity and V = I — VA^^(V-) is the Hodge projection 
operator to the space of divergence- free vector fields. For simplicity we assume 
that the forcing / is time-independent. 
In Fourier space, (|1.1[) can be written as 

(1.2) Vt + v\k\^v = -ikjPk [vjiv] + f , w(/c, 0) = vo, 

where v{k,t) = J- [v{-,t)] (k) is the Fourier transform of the velocity, i denotes 
Fourier convolution, a repeated index j indicates summation over j — 1,2,3 and 
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Pk = J-{V) is the Fourier space representation of the Hodge projection operator on 
the space of divergence-free vector fields, given exphcitly by 

(1.3) 

We assume that vq and / e and, without loss of generality, that the average 

velocity and force in the periodic box are zero, and hence v{0,t) = = /(O). 

Global existence of smooth solutions to the 3-D Navier-Stokes problem remains 
a formidable open mathematical problem, even for zero forcing, despite extensive 
research in this area. The problem is important not only in mathematics but it 
has wider impact, particularly if singular solutions exist. It is known that the 
singularities can only occur if Vw blows up. This means that near a potential blow- 
up time, the relevance of NS to model fluid flow becomes questionable, since the 
linear approximation in the constitutive stress-strain relationship, the assumption 
of incompressibility and even the continuum hypothesis implicit in derivation of NS 
become doubtful. In some physical problems (such as inviscid Burger's equation) 
the singularity of an idealized approximation is mollified by inclusion of regularizing 
effects. It may be expected that if 3-D NS solutions exhibited blow up, then actual 
fluid flow, on very small time and space scales, has to involve parameters other 
than those considered in NS. This could profoundly affect our understanding of 
small scale in turbulence. In fact, some 75 years back, Leray [22], [23], [24] was 
motivated to study weak solutions of 3-D NS, conjecturing that turbulence was 
related to blow-up of smooth solutions. 

The typical method used in the mathematical analysis of NS, and of more 
general PDEs, is the so-called energy method. For NS, the energy method in- 
volves a priori estimates on the Sobolev H™ norms of w. It is known that if 
t)||uii is bounded, then so are all the higher order energy norms t)||H'" 
if they are bounded initially. The condition on v has been further weakened [4] 
to /g II V X v{-,t)\\L'^dt < oo. Prodi [28 and Serrin |29j have found a family of 
other controlling norms for classical solutions I 21j . In particular, no singularity is 
possible if ||w(-,t)||Loo is bounded. The norm is also controlling, as has been 
recently shown in [31]. For classical solutions, global existence proofs exist only for 
small initial data and forcing or for large viscosity (i.e. when the non-dimensional 
Reynolds number is small). On a sufficiently small initial interval the solution is 
classical and unique. Global weak solutions (possibly non-unique) are only known 
to exist [22] , ^3 , [21] in a space of functions for which Vu can blow-up on a small set 
in space-timg^^l However, when / = (no forcing), a time may be estimated 
in terms of the ||wo||hi beyond which any weak Leray solution becomes smooth 
again. Such an estimate, which also follows directly from Leray's observation on 
the cumulative dissipation being bounded, is worked out in the Appendix. 

Classical energy methods have so far failed to give global existence because of 
failure to obtain conservation laws involving any of the controlling norms [32] . 



1-D Hausdorff measure of the set of blow-up points in space-time is known to be zero 

m 
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Numerical solutions to (II. ip are physically revealing but do not shed enough 
light into the existence issue. Indeed, the numerical errors in Galerkin/finite- 
difference/finite-element approximations depend on derivatives of v that are not 
known to exist a priori beyond an initial time interval. 

This paper introduces a new method in approaching these issues. In our formu- 
lation, the velocity v{x, t) is obtained as a Laplace transform: 

/•CO 

(1.4) v{x,t) ^ vo{x) + U(x,q)e^'^^^"dq, n>l 



Jo 

where U satisfies an integral equation (IE) which always has a unique acceptable 
smooth solution. Looking for v in this form is motivated by our earlier work |llj 
showing that small t formal series solutions, which exist for analytic initial condi- 
tions and forcing, arc Borel summable to actual solutions. In that case, the actual 
solution is indeed in the form (|1.4p . with n = 1. However, the representation (|1.4p . 
the IE for U{k,q) = !F[U{-,q)] (k). 1^1.6^ . and its properties important to (|1.4p are 
valid even when /, vq and the corresponding solutions are not analytic in x. An 
overview of our approach and nature of results is given in [12j . 

Note 1.1. For general initial data and forcing, U is in Li (M^, e^°"^dq), as defined 
in \2.8\} . If n > 1, then U is analytic in q in an open sector. For n ~ 1, the solution 
is q-analytic in a neighborhood ofM.~^ U {0} [H3 iff v{x,Q) and f{x) are analytic in 

X. 

As it will be seen later, using n > 1 is advantageous for some initial data. 
In Fourier space, (|1.4p implies 



(1.5) v{k,t)^vo{k)+ / U{k,q)e-'^^'"dq. 





Notation. Variables in the Fourier domain are marked with a hat Laplace con- 
volution is denoted by *, Fourier-convolution by J, while % denotes Fourier followed 
by Laplace convolution (their order is unimportant). 
As seen in §4, U satisfies the following IE: 

(1.6) U{k,q) = -tkj r g[q,q'-k)Hj{k,q')dq' + U'^''\k,q)=:N\lj\ {k,q), 

Jo 

where 

(1.7) Hjik, q) = Pfc [«oj*?7 + U.jivo + Ujlu] (k, q). 

The kernel Q, the inhomogeneous term U^'^\k, q) and their essential properties are 
given in (|422| and ([424|) in §4. 



Note 1.2. The solutions of ()1.6p . needed on M+, are very regular, see Note 11.11 
The existence time of v is determined by the behavior of U for large q. In this 
formulation, global existence of v is equivalent to subexponential behavior of U . 

The IE formulation was first introduced in pTj in a narrower context, and pro- 
vides a new approach towards solving IVPs. 



(3)This together with the estimate proves Borel summabiUty of the small t series. 
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2. Main results 



We define 



(2.8) 



Li(] 



dq) = < 9 



e-"«|g(g)Mg<oo 



Assumptions 1. In the following, unless otherwise specified, we assume that vq 
and f are in l^ i'o(O) = = /(O), n > 1, v > and a in h2. ^) is large enough 

(see Provosition 1 5. 1 /)) . 



Theorem 2.1. (i) Eq. kl.b]) has a unique solution U{-,q) E Li(R+,e °"^dq). For 
n > 1 this solution is analytic in an open sector, cf. Note [7771 We let U{x,q) = 

T-^ [t/(-,(7)] {x). 

(ii) With this U, v in \1.5\) (v[x,t) in respectively) is a classical solution 

of U.2\) ( U.l\) . resp. ) forte (0,a^i/"). 

(Hi) Conversely, any classical solution of v{x,t), t G (0,ro) has a Laplace 

representation of the form (77^ with U as in (i) and with 



U{k,q) [T[v{;T-'/-)][k)~T[vo][k)\ (q) 

a solution of HTg)) m Li(R+, e-"«dq), a > Tq"". 
The proof is given at the end of §5. 

Remark 2.1. Proposition 15.111 below provides (relatively rough) estimates on a. 
Theorem 13.11 gives sharper bounds in terms of the values of t/ on a finite interval 
[0, go]- Smaller bounds on a entail smooth solutions of (jl.ip over a longer time. 

We have the following result which, in a sense, is a converse of Theorem 112. 



Theorem 2.2. For f — 0, if lll.l]) has a global classical solution, then for all 
sufficiently large n, U{x, q) — 0{e^'-^"'^ ^ ') as q ^ +oo, for some C„ > 0. 

The proof is given in §7. 

Corollary 2.2. Theorems \2.1\ and \2.2\ implu that global existence is equivalent to 
an asymptotic problem: v exists for all time iff U decays in q for some n G 

The existence interval (0,a^^/") guaranteed by Theorem 12. II is suboptimal. It 
does not take into account the fact that the initial data vq and forcing / are real 
valued. (Blow up of Navier-Stokes solution for complex initial conditions is known 
to occur [30]). Also, the estimate ignores possible cancellations in the integrals. 

In the following we address the issue of sharpening the estimates, in principle 
arbitrarily well, based on more detailed knowledge of the solution of the IE on 
an interval [0, qo]. This knowledge may come, for instance, from computer assisted 
estimates or from rigorous bounds based on optimal truncation of asymptotic series. 
If this information shows that the solution is sufficiently small for q near the right 
end of the interval, then a can be shown to be small. This in turn results in 
longer times of guaranteed existence, possibly global existence for / = if this time 
exceeds Tr, the time after which it is known that a weak solution becomes classical. 
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(3.9) 



3. Sharpening the estimates; rigorous numerical analysis 
Let U{k,q) be the solution of (|1.6[) . provided by Theorem 12. II Define 

'U{k,q) for Q e (0, go] C M+ 







otherwise 



<.min{g,2qo} 



g{q,q';k)H^;^{k,q')dq' + m'Hk,q), 



Hf{k,q) = Pk 



Wo J* 



(fc,g). 



Using (|3.9p we introduce the following functionals of U^^\k,q), vq and /: 
(3.10) 6 := ai/2+i/(2n) / e-"«||;7(^)(-, g)||,i dg, 



(3.11) 

(3.12) 
where 



Bi 



go 



'B^{q')dq' 



i3i=4sup{|fc|So(A;)}||z)o||,i, So(A;) = sup Uq-q'f'^-"^^''^g{q,q'-M], 
fceZ3 qo<q'<q ^ ' 

B^iq) = 4 sup {|A;|Bo(fc)}||t/('^)(-,9)llii, ^3 = 2 sup {|fc|Bo(fc)}. 



Theorem 3.1. T/ie exponential growth rate a of U is estimated in terms of the 
restriction of U to [0, qo] as follows. 



(3.13) 



If ai/2+i/(2") > £^ 



2v eb then 



||?7(-,g)|hie-"«dg<cx3. 



The proof of Theorem 13.11 is given in 



Remark 3.1. In §8.1) it is shown that for a given global classical solution to (jl.ip . 
in adapted variables, the quantity ei + 2\/eb is small for large qo- 



Remark 3.2. In the proof it is also seen that if \\ij'^''\-,q)\\ii is small enough in a 
sufficiently large subinterval [qd, qa], then the right side of p.l3p is small, implying 
a large existence time (O, a~^/") of a classical solution v. The guaranteed existence 
time is larger if go is larger. If for / = 0, the estimated a"^/" exceeds Tc, the time 
for Leray's weak solution to become classical again (see Appendix), then global 
existence of a classical solution v follows. 



Since the improved estimates in Theorem 13.11 rely on the values of ?7 on a suf- 
ficiently large initial interval, we analyze the properties of a discretized scheme for 
numerical computation of U with controlled errors. 
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Definition 3.3. We introduce the following norm on functions defined on a S-grid 
in q 

= sup mi^i/"(5i-i/"(l + m252)e-""*||W^(.,m,5)||/i. 

Theorem 3.2. Consider a discretized integral equation consistent with ( cf. 

definition \9.2\) based on Galerkin truncation to [—N, N]^ Fourier modes and uniform 
discretization in q, 



see 77p below. Then, the error U — Ug^^ at the points q — mS, satisfies 
(3.14) \\U{;mS)~U^g''\;mS)y 
< 



1 ^am(5 



form > nig G Z+. where rrisS —: q„i > is independent of 5. In |g. Te,n is the 
truncation error due to Galerkin projection Vn and is the truncation error due 
to the S -discretization in q for a given N. We have ||Te,a,||("'''), \\{I-Vn)U\\^°'-^^ ^ 
as iV —> oo for any S and || J^f^'H'-"'''' — > as (5 — > 0, uniformly in N. 

Remark 3.4. For small q, independent of d, an asymptotic expansion of U exists, 
and solving the equation numerically for g G [0, qm] can be avoided. For this reason 
we start with q = qm- 



4. Integral Equation (jl.6p and its properties 
We define ii through the decomposition 
(4.15) v{k,t) = vo{k) +u{k,t). 

Then, (fL^ implies 
(4.16) 

Ut + J^jfepM = —ikjPk [vQ.jiu + Ujiv() + Ujiu] + Vi{k) —: —ikjhj{k,t) + vi{k), 

where vi is given by (|5.59p . Using v{k, 0) — vq, we have u{k, 0) = and we obtain 
from 



(4.17) u{k,t)=-ikj e^''l'=l'(*~'')/ij(fc,s)ds + t)i(A:) 



1 _ (.-•^\kft 



We look for ii in the form of a Laplace transform 

/•OO 

(4.18) u{k,t)= U{k,q)e-'^^^^dq; n>l 



We apply the inverse Laplace transform of (|4.17p with respect to r = 1/t" (justified 
at the end of the proof of Lemma 14.61 with more details in the Appendix) to obtain 
p.6p . The inverse Laplace transform of / is given, as usual, by 

(4.19) iP)^7r- f{s)eP'ds, 

ZTTl J c—ioo 
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where c is chosen so that / is analytic and has suitable asymptotic decay for Re s > 
c. 

For n = 1 the kernel Q is given by, see [llj . 



(4.20) g{q,q'; fc) = (Ji (2)^1(2') - Ji(z')>l(^)) 

z 

where z = 2\k\y/i7q , z' = 2\k\y^ vq', {n — 1) 
Ji and Yi are Bessel functions of order 1, and 

(4.21) U'-°\k,q) 2vi{k)^^^ , where z = 2\k\^ 

For n> 2 the kernel has the form (derived in the Appendix, see (|12.92p ) 
(4.22) 



1 

27ri 



c+zoo 



exp 



c— lOO 



i/|A:pT-i/"(l-s) + (g-(7's"")7 
(1 - s-")i/(2")-i(i _ s7i/")-i/V'/'J^(M) rfs. 



where 



(4.23) 



7 = - i^lfcl V/"(l - S7'/")(1 - s-")i/", 

9 



1 

27ri 



and C is a contour starting at ooe and ending at cx)e" turning around the origin 
counterclockwise. The function t/(°)(fc,g) in (|1.6p is defined by 



(4.24) U^'Kk,q)^^-^C-^{l 
where 

(4.25) G(A) - 



exp 



2TTi 



4.1. Properties of i^, G, Q, U^°'> and the relation between the IE and NS. 

Lemma 4.1. The functions F , G in Ii4.23\ ) and 14-25^ are entire andG'{iJL) — F{ii). 

Furthermore F{0) — yjjj:^, G{0) = and, for n > 2, their asymptotic behavior for 

large positive fi is given by 

(4.26) 

/ I _ 

3z7r 



7r(n + 1) 



r7,2("+i)^2(,x+i)Ini<^ exp 



F{fi) ^ { -I 



27r(n + 1) 



„2(„+i) ^2(„+i) exp 



2(n+ 1) 



27r(n + 1) 



„2(„+i) ^2(„+i) exp 



2(n+ 1) 



2{n+l) 



if arg fi = 



ifarg^e (0,^^^) 



ifarg^e (-^^,0) 
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(4.27) 



G{^x) ^ I 



where 
(4.28) 



7r(n + 1) 



/i2("+i)Im <^ exp 



1 



27r(n + 1) 



712(1+1) ^2(„+i) exp 



2(n+l) 



2Tr(n + 1) 



1,2(11+1) // 2(n+l) 



2(71+1) 
— ITT 



^2(r.+i) exp 



2(71+1) 



if arg ^ = 



ifargA*e (0,i^7r) 



ifarg^iG (-^71,0) 



Proof. These results follow from standard steepest descent analysis and from the 
ordinary differential equation that F and G satisfy, see ! jl2.1.1l | 

Remark 4.2. We see that F{fi) and G{fi) are exponentially small for large fi when 



, that is, when arg q G ( — — — ^ — — 



^ / (n— 1)^ (n— l)7r 

arg /it ^ ^ , J, u,,^^ 0.3, ouigy ^ ^ 2 ' 2 

Definition 4.3. For 5 > Q and n > 2 we define the sector 

jn- 1)tt ^ ^ jn - l)7r _ ^ 



Ss -.^ <q : a,igq e 



Lemma 4.4. For n>2,q,q'e e*"^K+ C Ss, with < \q'\ < |q| < oo and k £ I? 

we have 

where C2 only depends on 5. For n — 1, the same inequality holds for q,q' G IR+ 
with < q' < q. 

Proof. The case n — I follows from the behavior of Ji and Yi , see [UjiB- For 71 > 2, 
it follows from Lemma [4+1 that |/i^/^F(/i)| is bounded, with a bound dependent on 
S. Below, C is a generic constant, possibly S and 71 dependent. From (I4.22p we get 



(4.29) \g{q,q';k)\< 



I/V2|fc||g|l-l/(2«) 

x(i-.s-")i/(2")~i(i„5y/")-i/2 

For s G (l,i(l +7-1/")] we have 
(4.30) 

< c 



i(l+7-i/") 



,i/2|^||^|T3T7p^(^i+^2j; where 7 = 



l-l/(2n) 



(,_i)i/(2«)-iQ_iy/«^ 



-1/2 



1 + _ 1)1-1/(2")] _ i)l/(2„)-l (\_ 1^1/. 



-1/2 

2 2 

< C(l - 7l/«)-l/2 \l + {s~ l)l/(2")-l' 



(^)ln that paper the viscosity u was scaled to 1. 
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and for s e [5(1 + 7^^/"), 7"^/"), 

(1 _ g-n)l/(2n)-l(i _ s^l/")-l/2 < c[l + {s - l)l/(2n)-ll (1 _ s^l/")-l/2 

< C(l - ^l/")l/(2")-l(l _ s-^l/»)-l/2_ 



Thus 

/i < C(l-7i/")-i/2 



^(1+7-^/"), 



l + (5_l)l/(2")-l 
< C7-l/"(l - 7l/")-l/2 [(1 _ ^ (1 _ ^l/n)l/(2«) 

< C7"i/"(1 -7l/")l/(2")-l/^ 



/2 < C(l ~ -^l/«)l/(2n)^l 



i(l+7-i/") 



(l_sy/«)-l/2rf5 



< C7"^/"(l - 7i/")i/(2")-i/2^ 



I 



Lemma 4.5. (i) For n > 2 and =/= q Cz Ss, we have for a > I, 



C2V 



n/(n+ 



+1) 



llfc?7(°)(-,g)||,i <ci|||fc|f)i||,i|qri+i/"exp 
where ci and ci depend on S and n. Thus, we have 

(4.31) ^"e-"l«l||?7(")(-,g)||,id|g| < ci||^i||;ia-i/"r . 

PFif/i Ci = 1 and q G M+, i/ie bound in ^.31\ l holds for n = 1 as well. For n > 2, 

noting that Vi{0) = = /(O), we have 

(4.32) 



\\U^"H;q)\\i-d\q\<CG 



Vl 



iy\k\^ 



<CG<!|l^o||ii(l + -||^)o||/i)+- 



/ 



\k\' 



where 

poo 

Cg = sup n s-^\G{se'''')\ds. 

(a) If moreover jfep+^uo, \k\-' f G l^ (j — 0,1), then 

sup|g|i-i/"(l + |gne-"l«l||fc^"C/(")(.,g)||zi 

< 2ci|||fcKwi||,i < 2ci [,,\\\k\^+^vo\\ii + 2\\\k\HJii\\\k\vo\\ii + 

where the sup is taken over if n = I and over Ss if n > 1. 

Proof. The resuh follows from (|4.24p and (|5.59p using the asymptotics of G, ci. 
(|i?77)l and the behavior G(/i) ~ C/2 near fi = 0. For n = 1, the bound (|4.3ip 
follows from the fact that |2z~^Ji(2)| < 1. | 

The following lemma proves that a suitable solution to the integral equation 
(|1.6p gives rise to a solution of NS. 
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Lemma 4.6. For any solution U of S1.6\} such that \\U{-,q)\\ii e Li(M+,e °"^dq), 
the Laplace transform 

/>oo 

9/t" 



v{k,t)^Voik)+ / U{k,q)e-"i''' dq 
Jo 

solves U.^) for t g (O, a"^/") . For n = 1, v{k, t) is analytic in t for Re^ > a. 

It will turn out, cf. Lemma \12.1\ in the appendix, that \k\'^v{-,t) G l^ for t G 
(0,0;""'^/"). Therefore, v{x, t) = \v{-,t)\ (x) IS the classical solution of ULl]} . 

Proof. From (|4.24p . we obtain 

) 

e-''''"U^°\k,q)dq = Vi{k) 



-qt- 



27ri 



e'i^drdq 



1 — e^"!'^!^* 



Furthermore , we may rewrite (|4.22p as 



(4.33) g{q,q';k) 



1 

27ri 



1 pc-i-ioo 



J c—ioo 



l/"{exp -i^\k\^T-^/"{l-s) + {q-q'/s")T dr} 



ds 



since the integral with respect to r is identically zero when s £ (O, {q'/q)^^") (the r 
contour can be pushed to +cxd), we can replace the lower limit in the outer integral 
in Km by (q'/g)!/". Note that \\Hj{-,q)\\ii e Li (e-"l9ld|g|) , since 

(4.34) ||i^*G|U < ||^^|U||G||„ 

(see[Tn] and also Lemma [5.51 below). Changing variable g'/s" q' and applying 
Fubini's theorem we get 
(4.35) 

{k,q's'')Q{q^q',s;k)dq'^ds 



-zfc, / Hj{k,q')g{q,q';k)dq' 
Jo 

where for q > we have 











r 


ikjHj 









(4.36) Qiq,s;k) 



1 

27ri 



exp 



-i^|fe|V-i/"(l - s) + gr)] r-i/"dr. 



Laplace transforming (|4.35p with respect to q, again by Fubini we have 
(4.37) ^ e~'^*'" I^J^ {-ikjHj^{k,q's")Q{q-q';s,k)s"dq'ds^ dq 

= —ikj / ds g{t, s; k)hj{k, st), 



where hj{k,t) = C Hj{k,-) {t g{t,s\k) — C[Q{-,s;k)]{t "). By assumption, 

\\U{-,q)\\ii e Li {m+,e-"idq) and VQ{k) G l\ From (HHD and (lilM]) it follows that 
_f/j is Laplace transformable in q and 

hj{k, t) — Pk {voj*u + UjivQ + Ujiu\ {k, t), 
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while 

g(t, s; k) = t exp [~v\k\'^t{l - s)] . 

This leads to 



u{k,t)^t e-'^\k\ t{i~s) {k,st)ds + vi{k) 



1 - e 



(k, T)dT + Vi{k) 



and thus 



Ut + i'\k\'^u = —ikjhj{k,t) + vi , with 0) = 0. 



Therefore, using expression (|4.16p for hj, we see that v{k,t) — u{k,t) + Vo{k) 
()1.2p . with v{k,0) = vo{k). Analyticity in t of this solution in region Rej > a 
follows from the representation ()1.5p . It is clear that |/cp{)(-, t), f G l^, ensures that 
JF^i ^)] (x) is a classical solution to p.ip . | 



5. Existence of a solution to (jl.6D 
First, we prove some preliminary lemmas. 

Lemma 5.1. By standard Fourier theory, if v,w £ {'^^) > then so is viw, and 
\\v*w\\ii < I 

Lemma 5.2. 

\\Pk[w,iv]\\ii<2\\w,Mvy. 
Proof. It is easily seen from the representation of Pk in p.3p that 

(5.38) \Pkg{k)\<2\g{k)\. 
The rest follows from Lemma (|5.ip . | 

Lemma 5.3. Let C2 = C2{S,n) be given by 

0-2 = 2 sup Z/l/2|/fc||g|l/2|g-5'|l/2-l/(2n)|^(^^^/.^)| fom > 2, 

<j,<j'ee*'*R+C5a,0<|(j'|<|(?| 

C2 = 2 sup I/l/2|;.|^l/2|^^^^^,.^)| for„ = l, 

g,g'eR+,0<g'<9 

Then, for n > 2, we have 

f r\i\ ( 

(5.39) W](.,g)||,i < ^1/2,^,1/2 7^ (|<z|-.)-i/^+i/(^"){||?7(-,.e'*)|bi 

*||?7(-,se'^)||ii+2||i;o||,i||;7(.,se**)|hi}d,s+||J7(°)(.,g)||;i, 

(5.40) |lAA[C/W](.,(7)-AA[;7[2l](.,g)||,, 
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For n=l, 115. 39\) and hold for q G R+ , i.e. when (j) = 0. 

Proof. From Lemma l5.2| we have, for any q 

{k,q)\\ii < 2||;7(-, 9)11/1 * \\U{-,q)\\ii, 

and similarly 

\\Pk {voj*U{-,q) + Uj{-,q)*vo} ll/i < 4||^o||;i||?7(-,g)||;i, 

and ((09)) follows. 

The second part of the lemma follows by noting that 

(5.41) f/Wjf/W - J/f Si/I^l = U^t (f/W - U^'^) + {ul'^ ~ C>f ) tU^'^- 
Applying Lemma [5T2l to (|5.4ip . we obtain 

\\Pk{u\'^tu^'K-,q)-ulHu^'H-,q)}y<nu^^^ 

+ 2\\U^'^{;q)y*\\Um{.^q)-Um(^;q)y, 

from which (|5.40p follows easily. | 

It is convenient to define a number of different g- norms, q E e*'^R+ U {0} C Ss- 

Definition 5.4. (i) For a > 0, n > 2, we let be the set of analytic functions 
in Ss with the norm 

(5.42) supMi-i/"(l + |gne-"l^l||/(.,g)||/i<cx3, 

qeSs 

while for n — 1, A^"''^ will denote the set of continuous functions on [0, oo) with 
norm \\ ■ ||'"'. 

(ii) Let a > 0, 71 > 2, (5 > 0. We define a Banach space Ai''^ of functions along 
the ray \q\e^'^ € Ss with the norm 



(5.43) / e-"l^l||/(.,|g|e^*)||,id|g|<cx3. 

We agree to omit the superscript <j) when (/) — (which is always the case if n ^ 1). 
Lemma 5.5. We have the following Banach algebra properties: 

(5.44) \\ft9\\T^<\\f\\T%\\T'', 

(5.45) ll/S <Mo||/||(")||g||("\ 
where 

Proof. In the following, we take u(s) = ||/(., se*"^)]!;! and v{s) = |l5(., se"^)||/i . For 
(|5.44p we note that for any L > 0, 

pL n\q\ 

(5.46) / e~"l«l / u{s)v{\q\ ~ s)ds d\q\ 



e-""e-"(l«l-")u(s)w(|g|-s)dsd|g| < / e-^'u{s)ds ( e-'^^v{T)dT. 



Jo 
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From (I5.42p . we note that 

^i(s)^;(|gh5)ds< |l/|lM|l.g|l(")e«l«l ' 







- + s2][l + (|g| - s) 



21 



Finally, 

ds 



Sl-l/"(|g| - + 52][1 + (|q| - 



21 



_ + s2][l + (|q| - S) 

22-1/" Aq\/2 24-1/" /-"^ ds 



+ |g|2/4) 7o sl-l/«[l + s2] - + |5|2) sl-l/»[l + s2] ' 



where we used sup -^^l^^pj ,^ =4. | 



Lemma 5.6. Let C2 he as in Lemma \5.3\ and a > 1. The operator J\f in il.6]) is 
well defined on: 

(i) Ai''^ , where it satisfies the following inequalities 
(5.47) 

(5.48) \mu^^^] - U[U^^\T'^ < C2i^"'/'r(^i-^a-i/(2") 

X {(iic/W|ir^+ii?^'"iir^) iic/w-c/t'i|ir^+2iii)oihiii;7w -jjpiiir^}. 

(ii) A^"\ where it satisfies the inequalities: 
(5.49) 

(5.50) ||AA[C/[il] - AA[;7[21]||(°) < CaCai^-i/'a-i/f^n) 

X {Mo (||C/W||(") + ||C/[21||(")) IIC/W - + 2||f;o|hi||C/W - C/l2l||(")} , 

where C3 is defined in i5.52\) and depends on n alone. 
Proof, (i) For any < L < 00 and m > we have 

»l9l \ 

(k|-s)-i/2+i/('"^u(se'*)rfs 

1^1-1/2(1^1 _ 5)-l/2+l/(2n)g-o(|,|-.)^|^|^ 

< r e-'^'uise^'t') I /%'-'/'+'/('"^(s' + s)-i/2e-"^'rfs'Ls. 





(i 












>e''^)e-' 












Jo 



L 
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Using ((09l) it follows that 

/ e-^\'^^\mU]{-,\q\e'^)yd\q\ 
Jo 

<C2:.-i/2r(^i-)a-i/(2") (^[||?7||r*]' + 2||i;o||,i||C/||r*) +||?7(0)||f* 
From (|5.40p . it now follows that 

||AA[f7W] -AA[f7Pl]||,ie-"l9ld|g| 



< 



(ii) We first note that 

|5|l/2-l/n / e-"(l«|-^)(|g|-s)-l/2+l/(2")5-l + l/"(l + s2)-lrf5 

1 







|g|l/(2„) / g-a|g|(l-t)^-l + l/„(^_^)-l/2+l/(2„)(^^^2|g|2)-l^^ 







I /•1/2 |-l+l/n/-| j-\-l/2+l/(2n) 







1/2 I 



(5.51) < |g|l/(2")g-a|,|/2 Z"'^' + _ i)-l/2+l/(2n)^^ 







2l l/ri|g|l/(2n) rl i/2+l/(2n) 
l+klV4 71/2 ^ ' 

The first term on the right of ((53T|) is bounded by „2i/2-3/(2n)|g|i/(2n)g-a|9|/2 
For the second term we separate two cases. Let first a\q\ < 1. It is then clear that 

|g|i/(2") r g-akl(i-t)(i_i)-i/2+l/(2„)^^ 

7l/2 

^ 2n 



Now, if a|g| > 1, we have 



/■I /•1/2 
|^|l/(2„) / g-a|,|(l-t)(^_^)-l/2+l/(2n)^^^ |^|l/(2n) / g-a|,|*^~l/2+l/(2„)^^ 

Jl/2 Jo 

<W''-'r(i + i;)lo|,ir-'-'-'<o-"'->r(l + i^). 
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Combining these results we get 

r-l 



1/(2,1} 



1/2 



»kl(l-*)(l _ i)-l/2+l/(2n)^^ < 



where 

Therefore, 
(5.52) 



Ci = max <J r I ^ + 



2n 



2nJ n + 1 



sup 

kl>o 



\q\'-'/-{l + \q\^)e-'^\%\-'/^ 



l9l 



„ ^)-l/2+l/(2„)5-l + l/n(l + ^2^-1^^ 



l/(2n)p-7 4- 4«,iT.'v2+l/(2")^-7 



sup7"'^"'"'e ' + 4 sup7" ' ^"'"'e 



where 

Co = n2i/2-i/" 

L7>0 7>0 

From (|5.39p and the definition of || • ||("\ it follows that 

||AA[f7]||(") <C2C3^^-i/'«-i/(2») Mo(||C/||("))' + 2||f;o||/i||C>| 
Inequality (|5.50p follows similarly. | 



(a) 



1?7(°)| 



(a) 



Lemma 5.7. The integral equation ll.6\} has a unique solution in: 
(i) the hall of radius 2||f/(o)||"''^ in A'^'-'^ , if a is large enough so that 



(5.53) C2^-i/2r(^i-)a-V(2") (4||^o|ki +4||;7(°)||r^) 



< 1. 



Here C2 is the same as in Lemma \5.,S\ and depends on S and n for n > 2. For n = 1 
we have (j) — 0. 

(a) the ball of radius 2||f7(")||(") in A'-"^ if a is large enough so that 



< 1, 



(5.54) C2C3i^-i/'a-i/(2«) (^4||i)o||;i +4Mo||[/(°)||(")) 

where C2 (defined in Lemma \5. 3\) and C3 (defined in I15.52\) ) depend on S and n for 
n>2. 

Proof. The estimates in Lemma 15.61 imply that Af maps a ball of size 2 1 ^ m 
A"'^ back to itself and that Af is contractive in that ball when a satisfies (|5.53p . 
In A^"\ the estimates of Lemma [5.61 imply that Af maps a ball of size 2| 

|;7(o)||(a) 

to itself and that Af is contractive in that ball when a satisfies (I5.54p . | 

Remark 5.8. // a satisfies both h5. and \5.54^ , then it follows from Lemma 
and the uniqueness of classical solution of \1.2\ that the solutions U in Ai''^ and 
are one and the same. 

Lemma 5.9. The q-derivatives ofU{k,q) in A^"^^ for q > are estimated by: 



(5.55) 



dq' 



:U{;q) 



< Crn\\vi\\li 



-1+1/n -r. 



^aq+uia 



l + q^ 

where u) — q/2 for g < 2, = 1 for g > 2. 
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Proof. For q < 2, we use Cauchy's integral formula on a circle of radius q/2 around 
q and Lemma 14.51 to bound U for j^l > 0, argg G [— (?i — l)f + S,{n— 1)^ — S] 
(we may pick for instance S = j to obtain specific values of constants here). For 
q > 2, the argument is similar, now on a circle of radius 1. | 

In the following we need bounds on We rewrite (|1.6p using the divergence- 

free condition (note that kU is a tensor of rank 2) as 

(5.56) 

kU{k,q) = -^ik rg{q,q';k)Pk {Ujt[kjU] + Vo,j*[k,U]} (k,q')dq' + U^''-^\k,p) 



-■N 



kU 



where U^°^^\k,p) -.^-ik g{q,q';k)Pk Uji[kjVo] {k, q')dq' + kU^"^ {k, j). 

Jo '- 

We now think of U in (I5.56|) as known; then J\f becomes linear in kU. 

Lemma 5.10. // |fc|'^?}o G and a is large enough so that |5.5^[ ) is satisfied, then 

\\\k\W^ < 4ci [j^m^voWii + 2\\\k\vo\\l + + |||fc|i}o||,i. 

Proof. From (j5.56p . we obtain 
|||fc|f7||(") = ||A^[|fc|;7]||(") 

Lemma which applies when a satisfies implies that ||f7||(") < 2||;7(°)||(") 

and thus 

|||A:|f7||(") < C2C3i/-i/2a-i/(2")|||fc|[/|l("){2Afo||;7(°)||(") + \\vq\\ii} + \\U^'^'^^\'^"^ 

< ^|||fc|;7||("' + ||c/(0'i)||("). 

Thus, 
|||A:|;7|1(") 

< 2||C/("^i)||(") < 2|l|/s|?7("'|l(") +4MoC2C'3Z/-i/2a-i/(2")|||A:|vo|l,i|lJ7(")||("). 
Lemma follows from (|5.54p and bounds on U'^^^ given in Lemma [4.51 | 



II f II I f I 

Proposition 5.11. Assume /(O) = = wo(0) and we define || ||j — X]feez3\{o} fsp- 

If for n> 2, a satisfies the condition: 

(5.57) 



-1/(2") 



M\vo\ 



||{io||ii(l + -||?)o||;i) +- 





f 










/I- 


} 



< 1, 



with constants C'2 and Cq defined in Lemmas \5.3\ and \4.5\ then the integral equation 
U.6]) has a unique solution in a ball of size 2||[/'^'^)||"''^ in A"''^ . If in addition 
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jfcpwo G , then for n > 1 and a = ai is such that 

(5.58) C2^-i/2r(^i-)a-^/(^"'|4||€o||,. +4cir(^i)a7^/"||t)i||,.| < 1, 
where 

(5.59) vi{k) = {-i^\k\Hn - ikjPk [voj*vo]) + f{k) 

with ci defined in Lemma \4-.5\ then the integral equation il.b]) has a unique solution 
m a ball of size 2||C/(")||5'^^n A"''''*'. 



Proof. The proof follows from Lemma [4751 since (|5.57p and (|5.58p imply (|5.53p . and 
thus Lemma |5 . 71 applies . | 

Proof of Theorem \2l] 

Proposition l5. 1 II gives a unique solution to (|1.6p in some small ball in the Banach 
space Ai''^ for sufficiently large a. From Lemma 14.61 we see that U generates via 
(|f .5p a solution v to (|f .2p for t G [0,a^^/"). Classical arguments (presented for 
completeness in Lemma [12. II in the Appendix), show that |A:p'0(-, t) € l^ and hence 
[v{-,t)] (x) is a smooth solution to (jl.ip for t £ (0,a~^^^) . Analyticity in t for 
Re^r > ct follows from the Laplace representation. For optimal analyticity region 
in t, we choose n ^ 1. 

It is well known that (jl.ip has locally a unique classical solution [33], [H], [9j. 
Thus, given VQ,f S l^, all solutions obtained via the integral equation coincide. 
Furthermore, v{k,t) — vq is inverse-Laplace transformable in and the inverse 
Laplace transform satisfies ()1.6p . Therefore, no restriction on the size of ball in 
spaces Ai''^, A^"'' is necessary for uniqueness of the solution of (|1.6p . 



Remark 5.12. The arguments in the proof of Theorem 12.11 show that ||-0(-, i)||;i < 
oo over an interval of time implies that the solution is classical. This is not a new 
result. Standard Fourier arguments show that, in this case, we have ||?;(-, t)||Loo < 
00, i.e. one of the Prodi-Serrin criteria for existence of classical solutions [28], [29] 
is satisfied. 



6. Error bounds in a Galerkin approximation involving [—N, N]^ 

Fourier modes 

Definition 6.1. We define the operator M^^^ (associated to N ) by 
(6.60) N'^^^\ij\{k,q) 



ikj / g{q,q';k)VNPk 



(fc,g')dg'+7'Art/(°)(fc,g), 



where Vn, the Galerkin projection to [-~N,N]^ Fourier modes, is given by 



VnU 



(fc, q) = C/(fc, q) for k G [-7V, Nf" 
Lemma 6.2. The integral equation 



VnU 



(fc, q) = otherwise. 
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has a unique solution in yt^^fj as well as in , if a satisfies the conditions ir. 
Theorem \2.1i 



Proof. The proof is very similar to that of Theorem 12.11 part 1 , noting that the 
Galerkin projection Vn does not increase norms and A/'^^'' and Af have similar 
properties. | 

Lemma 6.3. Assume that a is large enough so that 

(6.61) ^2^3^.-1/2^-1/(2") (4|l{)o||,i +4Mo||;7(")||(")) < ^ 

and that \k\^i!Q, \k\f G l^ . Define the Galerkin truncation error: 



(6.62) Te,n ^VnU -Af^^^ 



VnU 



U 



VnU 



ikj / g{q,q'; k)VNPk Vo,ji{I - Vn)U + {I - Vn)Uj*Vo 



+ (I - Vn)Uj%VnU + VnUjUI - Vn)U Vn)U,1{I - Vn)U {k, q')dq' 



Then, 



where 



||;7_[/W||(") < ||(/-7'Ar);7||(") + 2||TB,Ar||("), 



||(/-7'jv)C/||("^+2||T£ 



N 



|(") 



< 



1 

N 

X 



2ci [^i^im-'voWii +2\\\k\vo\\ti + \\\k\f\\ii) + \\\k\vo\\ii 

{l + 4c||«o||;i + 12c (u\\\k\^voy + 2|| |A:|{;o||ii ||i)o||,i + ||/||ii)} . 

Proof. Clearly, 

||f7_[/W||(") < \\{I - VnWW^''^ + \\VnU - U'^^^l^'^l 
By (leleTj) . (|6?62)1 and contractivity of TV^^', 

<i||7,^[/_(7W||(") + ||Ts_^||W, 

so 

||;7_[/W||(") < ||(/_p^)f7||M + 2||T£,w||(«). 
Now estimates similar to (|5.49p imply that 

||r£.Ar||(") < c\\{I - P^)C/||(") [2\\voy + 2||7'Ar;7||(") + 11(7 - VnW"''' 

<c\\{I-Vn)U\\^"^ [2||i;o||,i+6||C/(")||(") 

and Lemma 15.101 implies that 
||(/_7,^);7||M<i.||fc,7||(") 
1 r 



< 



N 



2ci (i.|||A;|^i)o||,i +2|||fc|{;o||?i + +|||A;|i;o||,i 



(5)Recall this means A"-''' with (/- = 
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Hence the lemma follows. | 

7. The exponential rate a and the singularities of v 

We have already established that at most subexponential growth of ||[/(-, g)]];! 
implies global existence of a classical solution to p.ip . 

We now look for a converse: suppose has a global solution, is it true that 
U{-,q) always is subexponential in ql The answer is no. For n — \, any com- 
plex singularity tg in the right-half complex i-plane of vix^t) produces exponential 
growth of U with rate Re(l/ts) (oscillatory with a frequency Im(l/is)). 

However, if / = 0, we will see that for any given global classical solution of (jl.ip . 
there is a c > so that for any tg we have | arg ts\ > c. This means that for 
sufficiently large n, the function w(x,t~^/") has no singularity in the right-half t 
plane. Then the inverse Laplace transform 

1 rc+ioo 

can be shown to decay for q near 

We now seek to find conditions for which there are no singularities of v[x, r~^/" ) 
in {r : Rer > 0, r ^ M+ U {0}}. 

Lemma 7.1. (Special case of ^) // / = and v{-,to) E (T3[0,27r]), then 
v{x, t) is analytic in x and t in the domain |Im Xj \ < cv\t — tgl, < |i — to| < C* for 
arg(<: — to) € f ] J where c and C are positive constants (C depends on \\vo\\h^ 
and V and bounded away from when \\vq\\jji is bounded). 

See page 71 of [E]. 

Lemma 7.2. (i) Assume kvQ, f & l^ and a is large enough so that i f 5. 53\) holds for 
n—1. The classical solution of JJ. j|j has no singularity in Rcj > a, a; € T^. 

(ii) Furthermore, for / = (no forcing), no singularity can exist for arg(t — 
Tea) G for any < S < ^ and any x G T'^. (T^a is estimated in terms of 

ll'^^ollffi, J^, OLTid S in Theorem \12.1\ in the Appendix using standard arguments.) 

Proof, (i) The assumption implies Vq € H^iT'^). Since it is well known (see for 
instance [TS], [33], [S], [Hj) that a classical solution to (|l.ip is unique, it follows 
that this solution equals the one given in Theorem 12.11 in the form (II. 4p . From 
standard properties of Laplace transforms this solution is analytic for Rcj > a, 
where a is given in Theorem 12. II 

(ii) We know that under these assumptions ||w(-, t)||/fi ^ as t oo. There 
is then a critical time Tc^a so that standard contraction mapping arguments show 
that v{-,t) is analytic for t — Tc a & ds, seen in Theorem 112.11 in the Appendix. 
I 

Corollary 7.3. If f ~ 0, for any vq there exists a c > so that any .singularity ts of 
the solution v of is either a positive real time singularity, or else \ arg is | > c. 

Proof. If there exists a classical solution on R+ then ||t;(-, t)||/fi is uniformly bounded 
and by the proof of Lemma [721 (h) there is a Tc^a (as given in Theorem 1 12. II in the 
Appendix) such that v{-, t) is analytic for arg(i — T^a) S [— f j f ] • Let now Mi = 
maX(g[Q ||z;(-, i)||^i . Then by Lemma 1 7. II for any g [0,Tc_a + e] there exists 
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a C2 = C2(Mi) such that v is analytic in the region \t — t'\ < C2, \ arg(t — t')\ < 7r/4. 
Thus V is analytic in (see Fig[T]) 

{t:\t- t'\ < C2, 1 arg(i - t')\ < ^,0 < t' < T,.a + e} U {* ' I " ^ J}' 

Thus, if tg is a singular point of v, then tan | arg(is)| > c where 



Figure 1. The region of analyticity of v. 



C2/V2 



C2 



Proof of Theorem [272] 

By definition, 

1 



C/(A:,g) 



27ri 



-1/nA 



Vo{k) 



where the Bromwich contour C lies to the right of all singularities of u{k, t ^/") in 
the complex r-plane. By Corollary 17.31 u{k,t) has no singularities in the sector 

St.,p := |t : I argi I < := tan"^ c|, 

so u (fc,r-i/") has no singularities in the sector 

St,<p ■= : \ argr| < n^j. 

Clearly, if nc/) e (§,7r), then u{k,T^^/^'^) is analytic in a sector of width between 
TT and 27r, and in particular the Bromwich contour C can be chosen to be the 
imaginary axis. With the suitable decay of u{k, r^^/") at r = 00: 



u{k,t)=v{k,t)~vo{k)^0{t) asi^O, 

u{k, T 



which means that 



) = 0(r-i/") 



as T — > CXD, 
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Jordan's lemma applies and C can be deformed to the edges of the sector Sr,,/), i.e. 



U{k,q) = 



2TTi 



ooe-"**" Jo 



«(fc,T-l/")-{)o(fc)' 



dT 









2m 


[f ^ 


JO J 



t)(fc,T-i/")e«^ dr 



(carefully note that the integral of 'Co(fc) over the contour is 0). Further, as shown 
in Theorem 1 12. II in the Appendix, there is a sector S-^ in the right-half i-plane (with 
< 5 < f ) so that 



-Ti'Rot as t ^ C50 in 5*5. 



So 



||?)(-,r-i/")||,i < Ce-3''R<=(^ "'^) as T->0 along e±""*(0,oo), 
and the boundedness of ||'0(-, r~^/")||;i for large |t| implies that 

||w(-,r-i/")||zi < C7e-^"«'°(""''") for all r e e±™'^(0, oo). 
It follows that 



Jo Jo 



^i/r ' ^ cos (f)+qr cos n4> 



and a standard application of the Laplace method (with the change of variable 
r = g-"/("+i)5) shows that 

\\U{-,q)U < Cie-^=«''*"^" as q ^ +oo. 

8. Estimates of a based the solution of (II. 6|) in [0, go] 

Define J/^") as in (j?^ and U -U^°'\ Using (fTTB)) . it is convenient to write 

an integral equation for JJ^^^ for q > go- 



(8.63) 
where 



U^''\k,q)^-ik, / g(g,g';fc)ijf)(fc,g')dg' + f>(«)(fc,g). 



90 



min{(j,2go} 



t;(g,g';fc)i?j'^)(fc,g')dg' + [/(°Hfc,9), 



(8.64) U'-'\k,q)^-it 
and 

(8.65) i?j")(fc,g) -P,[z;o,,*m'') + i7]"'««o + t>j°^S?7("^](A:,9), 
(8.66) 

i7j'')(fc, g) = Pfe [i)o,,*?7(') + U^/^ivo + uftU^'^^ + ?7j'^S;7('^) + ufHu^'^] {k, q). 

Also, we define R^^\k,q) = —ikjH^-'\k,q). It is to be noted that the support of 
i?*^"-* is [0, 2go]. Thus, if J/^"-* is known (computationally or otherwise), then if^"^ 
and therefore U'^^^ are known for all g. 
Proof of Theorem [SHJ Note that 

|i?('')(fc,g)| < 2\k\ f2|{;o|i|?7('')| + 2\m-^\%\U^''^\ + \U^'^ \%\U^'^ ^ {k,q), 
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where | • | is the usual Eudidean norm in M'^. By Lemma [4.41 we can define a best 
constant 

(8.67) i?o(fc)= sup Uq-q'ff^-^'^^''^g{q,q'-M\ 

<?o<«'<9 ^ ' 

and conclude that 

\g{q,q'-k)R^''\k,q')\<2\k\Bo{k){q~q'Y/^^^^-''^[2\vo^^^^ 

+ \ij^'M\lJ^'^\] {k,q'). 

It follows from Lemma [5JJ that 



Biu + i?2 * w + B^u * u 



\\Q{q.q'-rW\;q')\\l^<H<l^q') 
where i}}{q) — g-^/^^")^^/^ and 

u{q) = ||C/(^)(-,g)|hi, = 4 sup {|fc|Bo(fc)}|| 



(<?'), 



i?2(g) = 4sup{|fc|i?o(fc)}||J7('^H-,g)lbi, 



B3 = 2sup {|fc|Bo(fc)}- 



Taking the /^-norm in k on both sides of (|8.63p . multiplying the equation by 
g-"9 fQj- some a > ao > and integrating over the interval [go, -M], we obtain 



M 



L,,M < / e-"« / ^J{q~q') 



Biu+B2*u+B3wu 



M 



{q')dq'dq+ "'^ u^'\q) dq 



< 



M , 



BlU + i?2 * u + * u 



M 



{q') / e-°'yj{q-q')dqdq' + / ""I u^'\q) dq 



go 

M 



< I e~°"^i^{q)dq / e-"« 
J qo 



Biu+B2*u+B:i,u*u 



M 



{q')dq'+ e-^'^u^'\q)dq, 



where 
(8.68) 



L 



10, M 



M 



e-'^'iu{q)dq, u(^\q)^\\U(^\-,q)\\i^. 



10 



If we use the fact that 



J go 



-aq 



u * v{q') dq' 



M 



qo 



u{s)v{q' — s) ds dq' 



qa 



AI 



qo 



AI 



u{s) / e'°"^ v{q' - s) dq' ds 

J s 

for any function v on [0, M] (recall that m = on [0, go]), then 



(8.69) Lg„,M < / e-"V(9)rf'Z 



Bi 



e-"« B2iq')dq 



+ &a-l/2-l/(2n)^ 
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6 = / e-"'^u^'\q)dq, 



90 



where 
(8.70) 

(8.71) ei = r 
For 

ei < al/2+l/(2«) _ ^l/2+l/(2n))2 ^ ^^^^ 

this leads to an estimate for Lq^,M independent of ill: 
1 " 

























Jo 













2ny 



.72) 



2e 



^i/2+i/(2„) _ _ ^(ei-ai/2+i/(2n))2_4e6 



So ||[/(-,(7)||ii e L^{e-°"idq) and the solution to ([ri|) exists for i e (0,a-i/"), if a 
is sufficiently large so that 

a > ao, ai/2+i/(2«) > ^ ^VTb. 

Alternatively, we may choose ao = a, in which case a has to be large enough to 
satisfy: 

^l/2+l/(2„) ^^^^ 

This completes the proof of Theorem 13.11 

8.1. Further estimates on ei, b and e. By Lemma 

(8.73) sup {|fc|9l/2(g_g')l/2-l/(2„)|t;(q^^'.fc)|| 
go<q'<q 

and by (|8.68p . I|8.64p . Lemma [4.41 Lemma H31 and the compact support of H'-°-\ 

(8.74) cs = sup |(7i/2-i/(2")y(«)(^)| < 

It follows that 



(8.75) 
where 



1 1 



b<CsT[- + —,aoqo], e < 2T[ - + — ]cgqQ ' , 



1 1 



2 2n 



-1/2 



r(a, x) 



is the incomplete Gamma function, and condition p.l3p is satisfied if 
(8.76) 



l/2+l/(2ri) 



-,1/2 



9o 



1/4 



If on a large subinterval [gc;,go]i ?)) ^^i^ therefore /f'"), decays, cf. the 

exponential decay in Theorem l2.21 then the estimated is small. Also, ei in (|8.7ip 
is small for large qo, ultimately since Bo{k) in (|8.67p is small. It is then clear that 
a in (|8.76p can be chosen small as well. 
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9. Control of numerical errors in [0, go] in a discretized scheme 

The errors in a numerical discretization scheme for 3-D Navier-Stokes cannot 
be readily controlled since these depend on derivatives of the classical solution; 
and these are not known to exist beyond some initial time interval. In contrast to 
physical space approaches, the q derivatives of the solution U to (|1.6p . are a priori 
bounded on any interval [qm, ?o] C for q„i > 0, by Lemma [5791 Further, if q„i is 
chosen appropriately small, then the small t expansion of NS solution provides an 
accurate representation for U and therefore of Hj in [0, qm] to any desired accuracy. 
Calculating the numerical solution to (|1.6p with rigorous error control is relevant 
in more than one way. 

In ^ we have shown that control of U on a finite q-interval provides sharper 
estimates on the exponent a, and therefore an improved classical existence time 
estimate for v. If this estimate exceeds Tc, the time beyond which Leray's weak 
solution becomes classical again (see the Appendix for a bound on Tc) then, of 
course, global existence of v follows. 

Furthermore, a numerical scheme to calculate (|1.6p . which is analyzed in this 
section is interesting in its own right. It provides, through Laplace transform, an 
alternative calculation method for Navier-Stokes dynamics. Evidently, this method 
is not numerically efficient to determine u(x, t) for fixed time t; nonetheless it may be 
advantageous in finding long time averages involving v and Vu needed for turbulent 
flow. These can sometimes be expressed as functionals of U. 

Definition 9.1. We introduce a discrete operator Afg^^ by 

m — l 

(9.77) {aA]^^[T>]} (fc,m(5) = w^^\m,fn';k,S)VNH^J\k,m'6) 

+ Lf^°'^\k, mS) - ikjw'^^^^^ (to, k, 5)VNH^^\k, m6), 
where k G [— iV, N]^ \ {0}, N 9 m > TOs, q-m = (qm is independent of S) and 

(9.78) U^°'^^ (fc, mS) U^"^ (k, m5) - ^% Z"^"' ^(m<5, g'; k)VNHf'^ (k, q')dq' 
is considered known, while for to' > rris, 







(9.79) 

H!j^g\k,m'S) = Pk \vo,j{k')V{k - k',m'S) + vo{k')Vj{k - k',m'5) 



k'e[-N,NY\{QM} 



+ Y Pk\yj{k',fn''5)V{k-k',{rri' -rn'')5) w^'^\rn',rn''-,k,5) 



m —ms..,m —rris 



k'e[-N,Nf\{0.k} 

rris — 1 

+ 2 Y w^'^^^\m' ,k,5)PkW^\k)iV{k,{m' -1)5) 



1=0 



In (|9.79p . E^^\k) involves 'Do(fc)-this representation encapsulates the singular con- 
tribution of U{-,q') and U{-,q — q') when q' and q — q' are small respectively. The 
precise form of these functions and of the weights w^^-* (to, to'; k, S), w^^-'^\m,k,S), 
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y;(2) (^jjj^' ^ a^j^fj ^(2'') (m', fc, (5) generally depend on the particular discretiza- 

tion scheme employed to calculate A/'j^'' [U]. Also, note that in (|9.79p , the nonlinear 
terms in the summation are absent when mg < m' < 2™^. To simplify the discus- 
sion, we do not specify the weights, but only require that they ensure consistency, 
namely that in the formal limit 6 0, the discrete operator A/j^^ becomes Af^^\ 
Based on behavior of the kernel Q, consistency implies that 

C Ai/(2") 

(9.80) |fc||w;(i)(m,m';fc,<5)| < ' 



mi/2(TO-r7i')i/2-i/(2") 

\kw^^'''>\ < Ci^i(5l/2+l/(2n)(^J)-l/2 ^ 1^(2) I < ^^g^ 1^(2,0 I < C^S^/" (1 + 1)-^ + ^" , 

Consider the solution 

(9.81) U^'^\k,md) = jw]^^ U^^^ I {k,mS) for < m, fee [-iV, Af]^, 

where as noted before, qm = ni-sS is small enough so that the known asymptotic 
series of t/ at g = can be used to accurately calculate t/'-^' and -ffj^"* for q < qm, 
and thus of ;7(°'^) and in (p78)l and (19791) . 

Definition 9.2. We let 

&e t/ie truncation error due to q- discretization for a fixed number of Fourier modes, 
[—N,N]'^. The discretization is consistent (in the numerical analysis sense) ifT^^g 
scales with some positive power of S and involves a finite number of derivatives of 
U. 

Definition 9.3. We define \\ ■ H^"'''-', the discrete analog of \\ ■ \\^"\ as follows: 
sup mi-i/"<5i-i/"(l + m2j2)e-"™5||/(.,TO^)||ji. 

Remark 9.4. More specific bounds on the truncation error depend on the specific 
numerical scheme. It is however standard for numerical quadratures to choose the 
weights w*--'^ so that ^-integration is exact on g S [g'mi'Zo] for a polynomial of some 
order I. For a general V{-,q), the interpolation errors involve / -I- 1 q-derivatives. 
Lemma lS. 91 guarantees that the derivatives of U are exponentially bounded for large 
q. It follows that WT^^^sW^'^'^^ as (5 ~> 0. 

Remark 9.5. In the rest of this section, with slight abuse of notation, we write * 
for the discrete summation convolution in q-space ( i. e. sum over m' ) and * for the 
discrete double, Fourier- Laplace, convolution. Since the rest of the paper deals with 
discrete systems, this should not cause confusion. 

Lemma 9.6. For m > rus, Hj^\-,m6) satisfies the following estimate: 

(9.82) \\H^J\;m5)y 



< C- 



for some known constant Ce 
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Proof. Using the properties of discrete convolution we see that 

<chv,\\i4uf\-,ra5)\\i^+5 ^ ||t/f' (•, m'(5)||,i IjC/f (m - m')<5)|hi 

m'—nis 

nis — 1 

+ S'/" (m' + l)-i+i/«p('"')|hi||C/f ^•,(m-m')<5)||zi 
where, by a standard integral estimate, 

m — 1 

' m ' (1+m d I > ^ 5- < 

„5^i [5m'S{m - m')]i-i/"(l + (52m'2)(l + 6^{m - m'f) 

nis — l f. 
rl-l/n l-l/rw-, I 2 r2N ^ 

^^J5(m' + l),5(m-m')]i-i/"(l + (S2(m-m')2) 
for C independent of m, to' and (5. In the above estimates we have used 

which can be obtained from the definition of i?*^'" ^ . | 

Define H^-^'^'' and i?-^'^'' by substituting f/j^'' = t/]^'"'^'' and Ug^'^\ respectively, 

inHj;). ' 

Lemma 9.7. For m > m^, we /laue 



- (i + to2,52)to1-i/«5i-i/"||''^ II 

x{||;7f'^)||(">) + ||?7f-^)||(">) + ||«oll..+c^}. 

Proof. The proof is similar to that of Lemma 19.61 | 

Lemma 9.8. (i) For C4 defined in i9.85\) . assume a is large enough so that 
(9.83) 2C4a-i/2-i/(2«) ^ ||^^||^^] ^ 2||C/("^^)||("'*)) < 1. 

Then, for any (yT^ > (^o > > 0, A/'j^' is contractive and there is a unique solution 



, which satisfies the bounds 



('ii^ //a «s such that 

(9.84) 2C4a-i/2-i/(^") ([C^ + ll^olb^] + 2\\m''^Y"''^) < \, 
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then 

Proof, (i) Wc have 

(9.85) IIAAf )[;7n(-,m<5)||,i < ||C/(°'^)(., m5)||,i 

m'—rris 

< f J \\fj(0,N)u(a,5) 

- (l + m2,52)TOl-l/«51-l/n\" " 

+ C4a-V2-i/(2")rri|(-^)(||C/ri|(-^) + ||^)o|hi + 
where, by a standard integral estimate, 

' ^ [M^/^[<5(m-m')li/2-i/(2")[5TO'li-i/»(l + TO'2,52) 

m'—ms 

and 

^l/2+l/(2n) ,/2+l/(2n) < -l/2-l/(2n) 



(AT) 



has a unique solution such that 



Hence the first part of the lemma follows, 
(ii) Under the assumption, 

ll^(A) _ j;>W|j(o,5) < ||aaW[;7W] _ AAf)[f7f )]||(«.*) + |irW||("-^) 

So 

and the second part of the lemma follows. | 
Proof of Theorem [HI} Note that 

From Lemmas 16.31 and 19.81 it follows that 

which tends to zero as iV ^ oo, (5 — > 0, by Lemmas 16.31 and 19.81 
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10. Numerical Method 

In this section we describe a numerical scheme for calculating the solution C/j^'' 
over a fixed interval. The procedure can be further optimized in a number of ways, 
such as adapting the quadrature scheme to the features of the kernel. 

10.1. Outline of the Algorithm. The main algorithm is summarized as follows: 
initialization; 

startup routine; 
for each time step 

advance the solution using second order Runge-Kutta integration; 
end 

estimate the error and output the results. 

10.2. Startup Routine. One difficulty in numerically solving (|1.6|) is that the 
equation is singular at g = 0. To overcome it, we first compute ii for small t by 
solving (|4.16p using Taylor expansion: 



m— 1 



where 

ci = -Oi, 

v\k\^Cm - ikjPk ( V0,j*Cm + CmjJwo + ^ CljiCm-£ 



1 

Cm+1 



m — 1 



m + 1 

Then U is computed for small q by 

mo 



=1 



m > 1. 



where 



r(m/n) 

10.3. Second Order Runge-Kutta Integration. After computing the solution 
on [0, (7„i] for some q„i > by using Taylor expansions, we solve the integral equa- 
tion (|1.6p on [(7m, go] using second order Runge-Kutta (predictor-corrector) method. 
Since this numerical scheme is preliminary and far from being optimized, we do not 
include the details here. 

What is worth mentioning is the evaluation of the functions F(^) and G(/i). 
As shown in earlier sections, both F and G are entire functions and have power 
series expansions at /i = 0. For small /x, these expansions converge very rapidly 
(super-factorially) and provide an efficient way to evaluate F and G. For large /z, 
however, the alternating nature of the expansions raises the issue of catastrophic 
cancellation, and it is no longer appropriate to use them for numerical computation. 
In this regime we use the asymptotic expansions of F and G, which we derive below. 



that 

roo roc / m\ 

so T{m/n) ■ dm = Cm- 
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While the complete asymptotics of F and G can be derived using Laplace's 
method, a faster and easier way is to use the differential equations they satisfy. For 
example, recall that for n = 2, 



where 



Ii=i 



It is easy to check that h satisfies the third-order ODE (the same equation satisfied 
hyF) 

+ 1'{ - 2/1 = 0, 

and it has the leading order asymptotics 



7i ~ 2 



where 



= 3-2-2/V/V'^/^. 
If we make the change of dependent variable 



then Ji must have the form 



7i=2^/-e-Vi(^), 



m— 1 



and it solves the ODE 



It follows that 



jr-3J!' + (3 + ji^)ji-jLA=0. 



where ai, 02, etc. are determined by the recurrence 

ao = 1, 

(^12m^-12rn+ljaTO_i + (aiti^ - 1 2m^ + 9m - 2 j a„ 



ai 



1 

12' 



1 



12m 
Similarly, 



m > 2. 
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where 

Co = 1, 



Cl 



24m 



12' 



C2 



35 

'288' 



-48771^ + 60m - 2^c,„_i + (^-32m^ + 108m^ - 80to + 9^Cm-2 



(^-8m,^ + 52to^ - 102to^ + 67m - 14^ c„ 



m > 3. 



11. Preliminary Numerical Results 

For all computations in this section we take n — 2. The numerical results 
and computation scheme are preliminary. The algorithm has not been optimized 
for efficiency, and not all estimates have been rigorously analyzed yet, and these 
will be published elsewhere. Nonetheless, the partial results show some important 
features of the integral equation approach. 

11.1. Test Case. We first tested our code with the following test function: 



(Kida flow) : v = {v'-^\ v^^\ v^^^ ) , 



smxi 
l + t 



(cos 3a;2 cos X3 — cos X2 cos 30:3), 



V^^\xi,X2,X3,t) = v'^^\x3,Xi,X2,t) = v'^^\x2,X3,Xi,t). 



The forcing / corresponding to v was generated with v = 1 and equation (|1.6p was 
solved without the knowledge of v. The computed solution was then compared to 

V. 

For this test case, the startup routine computed the solution on [0, g™] = [0, 0.2] 
using mo = 8 terms and the Runge-Kutta solver advanced the solution to 90 = 1- 
2N = 16 points (i.e. 8 Fourier modes) were used in each dimension (excluding the 
extra points for anti-aliasing) . 

We computed the solution for different step size S and the errors at go 

eg = ma.x\U^g^\x,qo) - U{x,qo)\ 

are listed in Table [T] To ensure the error decays at the right order 0{S'^), we also 
included in the table the numerical order of convergence: 

f3s = log2 — . 



Table 1. Test case: errors at qo- 



5 


es 


I3s 


1/20 


1.3399e-04 




1/40 


3.1987e-05 


2.07 


1/80 


7.1462e-06 


2.16 


1/160 


1.3620e-06 


2.39 
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11.2. Kida Flow. Now we consider the Kida flow with the initial condition 

Vq^'' X2, a;3, 0) — sina;i(cos 3x2 cosxs — cos X2 cos 3a:;3). 

We computed the solution for j/ = 0.1 with zero forcing to go = 10 using 2N = 128 
points in each dimension, and step size S — 0.05. The parameters for the startup 
procedure are the same as before: Qm = 0.2 and mo = 8. To investigate the growth 
of the solution J/]^'' with q, we computed the ^-'^-norm 

kel-N,N]^ 

and plotted ||?7f ^^(•, vs. q in Fig[2] For comparison we also included in Fig[2] 
a plot of the solution to the original (unaccelerated) equation. 

FigE] shows the plot of log||C/f''Y,g)||ii vs. q^/^. Note that \\U{-,q)\\ii - 
cie-=2«''' for large q, where ca = (0.3)2/32-5/33 ^ o.42. 

11.3. Longer Time Existence. We next computed the constants in estimate 
(|8.69p . By taking go = 10 and oq = 30, we obtained 

6«0, 1.1403, ei« 13.6921. 

This implies the existence of the solution for a > 32.7564, which corresponds to an 
interval of existence (0,0;-^/^) = (0,0.1747). 

We compare with a classical estimate of the existence time. The formula 

T 1 

(where c„i is known) was optimized in the range m > 5/2, giving a maximal value 
Tci « 0.01 at m « 3.2, about 17 time shorter than the time obtained from the 
integral equation. 

Furthermore, considerable optimization of code is expected to allow numerical 
calculation over much larger q-interval. 

12. Appendix 

12.1. Derivation of the integral equation and of its properties. 

The integral equation. We start with the Fourier transformed equation (jl2.86[) : 

(12.86) ut + i'\k\'^u — —ikjPk[vo,j*u + Uj*vo + ujiu] + vi{k) 

= : —ikjhj + vi{k) =: f + Vi{k), 

u{k,0) = 0, 

where 

vi{k) = f{k) - i^lfcp-Co - ikjPk[vQj*vo]. 
For n > 1, look for a solution in the form 

POO 

(12.87) u{k,t)^ U{k.,q)e-'i/*" dq 

Jo 
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where 

(12.88) r{k,t)^-ikfhj{k,t)=-ikj I Hj{k,q)e-'^^^" dq 



POO 

= : / R{k,q)e-'i/'" dq. 
Jo 

Inversion of the left side of l|12.86p and the change of variable t ~ t~" yield 
(12.89) u{k,t)^ / e-''l''^l'(*-^)f(A:,s)ds+ / e~''^''^''^'~'Hiik) ds 



^ f te-l^l^'(i-)f(fc, ts) ds + ^ f 1 - e-I^I^O 
Jo ^ ' 

' r-i/"e--l'=l'^"^'"(i-^) r R{k, q')e-^'^-"- dq'ds+^^(l- e-H^I^--^'") 
Jo nk\ ^ ' 

= ■ I{k,T)+J{k,T). 

Inverse Laplace transform (formal for now) of / and J yield: 
(12.90) — / I{k,T)e'i^dT 

ZTTZ J c—ioc 

R{k,q') I— T-l/«e-HfcP--'^"(i--)+(W.-")-dr dsdq' 

Jo I Jc-ioo J 

/•OO 

= : / R{k,q')g{q,q';k)dq', 
Jo 

where 

C = (g - g's-")r, M = - s)(g - g's"")!/", 

while 



:y\k\^q\2TTt J,_, 
where 

The Bromwich contour is homotopic to a contour C from ooe^*'^ to the left of 
the origin, ending at ooe™ encircling the origin, and we finally obtain the integral 
equation: 

U{k,q) - r giq,q';k)\-ikjHj{k,q')\ dq' + U^''\k, q), 



where 



Hj{k,q) = Pk 



vo,i*U + Univo + U^IU 



(fc,g). 
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Rescaling the integration variable, s 57^/", the kernel in ()12.90p becomes 



(12.92) g{q,q';k)^q'/''-'-/'/" I {I - s-'^f'^-^ F{^i) ds 

(1 - s-")i/(2")-i(i _ s7i/")-i/V'/'J^(/")ds, 



i^i/2|fc|gi-i/(2«) 
where 



and 



Furthermore, from (|12.9ip we have 

(12.93) {7("'(fc,g) = ^G(Hfc|\^/"), 



where 



G{^l) = -^. f e'^-^i-''" dC, G(0) = 0. 



Power series representations of F and G. To show that F is entire, we start with 
the definition 

(12.94) F(m) = ^ / C"'/"e«e-^«"''" dC 

J c 

and expand e~^^ into power series of C^^/" to obtain 

where the interchange of order of summation and integration is justified by the 
absolute convergence of the series. From the integral representation of the Gamma 
function (see [1]) we get 



J ^-0+i)/«eCdC = 2isin(^^7r^r(^l- 



r((j + i)/n)' 



(where in the last step we have used the identity sin(7rz)r(l — z)T{z) — tt) and thus 
F has the power series representation 

°° f— IV 

F(/x) = ^Fj/Lt^ where Fj - ^ ' 



.^0 mij + i)/ny 

Similarly, G is an entire function and has the power series representation 

(-1)^- 



G{fi) — ''^Gjfi^ , where Gj 
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12.1.1. The Asymptotics of F and G for n > 2 and large jj, > 0. Elementary 
contour deformation and estimates at show that 



where 

(12.95) h{li)= r-i/"e''^/" exp 
Jo 



-r - /ir-i/"e*''/" 



dr 



= nMi-2/("+i)e''^/' 



Jo 



s" ^exp 



_^»/("+i)(s" + e'''/"s-^) 



ds 



where 



Similarly, 



w 



/•OC 

= „^i-2/("+i)e2i'^/(»+i) / a;"-" exp -w<f{x) 
Jo L 

= ^"/(»+i)eW("+i)^ ^(a;) = a;" + i. 



da;, 



/■OO 

i7r/(n+l) / ^n- 

Ja 



' exp 



—wip{x) 



dx. 



We now use the Laplace method to obtain the complete asymptotic expansion of 
Ii for large w with argw € (— f , f ) or arg/x e ^— ''"2^''^ , ^"2to"'^ ) ■ then show 
that /i solves a linear differential equation. It will follow, from standard results on 
asymptotics in ODEs, that the expansion is valid in a wider complex sector. First, 
it is easily seen that the only solution to the equation 

if'ix) = na;"-^ -4=0 

on the positive real axis is a; = a;o = n~^/^"+^^. If we introduce a new variable 
(12.96) C = 

then clearly ^ decreases monotonically from a; = 0+ to a; = a;o, where it attains the 
minimum value 

eo = ^(.xo)=n-"/("+i)(n + l). 

We denote this branch of ip^^ by a;i(^). Further, as x increases beyond a; = a;o up 
to 00, ^ increases from ^0 to 00. We denote this branch of ip~^ by a;2(^). It follows 
that 



I, = „^l-2/(n+l)g2W(n+l) 



To find an expansion of a;i(^), i = 1, 2, we note that 



d^+ 



d^ 



J=2 



and thus 



00 ( _ \j 1 

(g-Co)-E^''^(^o) / = i^f"{x,){x-x,f, 

j=3 
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or 



(12.97) 



X± — Xq ± 



3=3 



where X- = Xi and 2;+ = X2- By (|12.97p we have 



nxr\o-^7'\o 



Watson's lemma then implies that 
(12.98) 



„^l-2/(n+l)g2i7r/(«+l)g-5o«' 



OG / 



We see that 



/,[2] _ .[1] _ I even 
^ ^ ~\2bf J odd 



Similar analysis for Ii gives 



(12.99) Ji ~ V-2/("+i)e"2-/(n+i)g-Co^ ^ ^ (6f -feW)V/2e-'="rf77 



i-2/(n+i)g-2W("+i)e-?oS' ^ (^5(2] _ ^[i]^ r( 1 + ^ I ii; 



p ,T,-l-j/2 



With ^o^'^ replaced by z, we finally obtain for fj, large and positive 

(12.100) F{fi) = -Ini/i 

TT 

~ ^Ini|;,(»-2)/[2(n+l)]g3W[2(n+l)]g-. ^ 26[l_ir(^m+ ^ 

where 

(12.101) ^0 = n-"/("+i)(n+l), z = ^oM"/("+i)e*'^/("+i). 
A similar analysis shows that 

where z, ^0 are given by (|12.10ip and d'*' are coefficients of the expansion 

xrHo 
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To obtain the leading asymptotics of F and G, we note that 

n-2 n-2 n-2 



^ - x^"^ </' 
It follows that 

V293"(xo) V2 
Similarly 



n— 1 -1 



As a result, we have to the leading order, 



G(;i) ^ -y|nl/[2(»+l)l(„ + l)-l/2lni|//V[2("+l)]gW[2(n+l)]g-.y 

Differential equations for F and G for n € N and extended asymptotics. To derive 
a differential equation satisfied by F, we differentiate (|I2.94p n times in fi (justified 
by dominated convergence) 

Integrating by parts once, we get 

2Tn Jc \ n J 

= i-irnFif,) - f,F^-+'\f,), 

so the differential equation satisfied by F is 

Since G' — F, the differential equation satisfied by G is 

/iG("+2) + G("+i' - (-l)"nG' = 0. 
Integrating once and using G(0) = 0, we obtain 

(12.102) /iG("+i) - (-l)"nG = 0. 
We can make the same argument for 

.0+ 

(12.103) h{ti)= / e'^-'"^'""dC-l 



or for 

(12.104) h{li)= e^-"^ ''"dC-1. 



OQC^ 



It is to be noted that G{fi) = — , while l2{p) = h{^) and /2(a^) = 

hiti). 
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Equation (I12.102P has (n + 1) independent solutions with the following asymp- 
totic behavior for large /i (see |34]): 
(12.105) 

' ; z ^oe*'^/("+i V"^^"'^^^ , j = 0, 1, . . . , n. 



Thus, there is only one solution with the asymptotic behavior 

^„i/[2(n+i)](„ + i)-i/2^«/[2(n+i)]g^p[_^] for arg z = 



(all solutions independent from it are larger). Since I2i^J■) has this asymptotics in 
particular for arg/i = — ^, corresponding to argz = as discussed already, I2 is 
the only solution of (I12.102[) satisfying 

(12.106) /2(A*)-^-\/^ni/[2("+i)l(n + l)-i/V"^f'*"+'^'cxp[-2] for arg// = --. 

V TT n 

As we rotate around in the counter-clockwise direction starting from arg z = in 
the complex z (or complex fi) plane, the classical asymptotics of I2 can only change 
at antistokes lines. The first antistokes line is argz = f + ;^fri corresponding to 

(n+3)7r 

Similarly, in a clockwise direction, the first antistokes line is arg z = ~ , 

I.e. arg/i = - ' ■ 

Therefore, for arg/i G ^'"2,^"'^ ) ''"2rt"'^ ) ^-'^'^ asymptotic expansion /2 is the 
same. 

From the symmetry between I2 and I2, it follows that 

(12.107) /2(/.)~-y|ni/P("+i)l(n + l)-i/V/[2(«+i)lexp[-z] 

for arg/i G (-^^, ^^). Since G(//) - ^i- [hifi) - /2(/^)] , noting that hifi) 

is dominant for arg/i G ^0, , it follows that in this range of arg/i, G(/i) '--^ 

— 27-^2 (/i)- while for arg/i e (^ ^"211^^ '^)' ^i^^e /2 is dominant, G{p) ~ ^hifJ-)- 
Lemma 14.11 follows. 



12.2. Instantaneous smoothing. The following result shows that the solution 
v{k,t) obtained from U{k,q) corresponds to a classical solution of (|l.ip for t G (0, T] , 
i.e. there is instantaneous smoothing due to viscous effects. This is a known result 
(See for instance |5]), but we include it for completeness. 

Lemma 12.1. Assume wo,/ G where va{0) = = /(O). Assume further 

that il.l]) has a solution v{k,t) with ||u(-,t)||;i < 00 fort G [0,T]. Then v{x,t) = 
!F~^ [{)(•, t)] [x) is a classical solution of for t G (0, T]. 

Proof. It suffices to show |fcpi)(-, t) G l^ for t G (0, T] since this implies w G and 
usual arguments imply that v satisfies (jl.ip . 

Consider the time interval [e,T] for e > 0, T < a^^/". Define 

'We{k) = sup \v\{k, t). 

e<t<T 
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Since \v{k,t)\ < \U{k, q)\e^°"^dq, wq (or We) satisfies 

/"OO 

\\woy < / ||C/(-,g)||zie-"«dg. 



On [e,T] for e > 0, (HH]) implies 
(12.108) 

vik, t) = ~ik, [ e-''l'=l'(*-^)Pfe (vjiv) (fc, r)dT + «oe-'^l'=l'* + ^ (l - e-^'I'^'l'*) . 

Jo '^I'^l ^ ' 

Therefore, 

\k\\v\(k, t) < 2 {woiwo} f |fc|2e-''l'=l'(*-"'dT + |fc|{)oe-''l*'-|'* + ^17(1- e"''!*'-!'*) . 

Jq v\k\ \ I 

Since i/j/cpe-'^l'^l'^'-^^dT < 1, it follows that 

(12.109) |fc|w,/2 < - {wo*i&o} + \/ — ( sup7e-^' ) |{)o| + ^ 

V J^e V7>o / v\k\ 

Using now the bounds on Wq we get 

/ 



e/2lhi < ^ {r(-,g)ll?}' + ^TT^II^ollz. + 



lfcl 



The evolution of v is autonomous in time, and thus, for to. [§, P] we have 
(12.110) t)(fc,<) = / e-''l''l'(*-^)Pfe {vf*{kjv\) {k,T)dT 

Je/2 

+ ^(fc, e/2)e-H'=l=(*-/2) ^ ^ 

where we used the divergence condition k ■v{k,t) = 0. Multiplying p2.110p by |/cp 
and using (|12.109p . it follows that for t E [e,T] we have 

\k\^\v{k,t)\<2w,/,i[\k\w,/,] f |A;pe-H'=l'(*--)dT 

Jt/2 



(sup7e-'' ) |v(/c,e/2)| + 



t/(t - e/2) V7>o y /I ' ' 

implying that 

IllfcP^i.llii < l\\w,/2\\i4\k\w./2\\i- + §ll^./2lhi + 
Since e > is arbitrary, it follows that \k\'^v(-,t) e for t e (0,T]. | 

12.3. Estimate of Tc beyond which Leray's weak solution becomes classi- 
cal. It is known that (|1.2p is equivalent to the integral equation 

rt 

(12.111) v{k,t)^ e-''l'=l'(*-^)Pfe [-ifcjf)j«{;] (fc,T)(iT + e-''l'=l'*f)o 

^:F{M[v]{;t)}ik). 
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Applying T^^ in fc to ()12.11ip . it follows that 



(12.112) v(x,t) = e^^^VQ~ \ e^'^'-^^^V{{v-^)v\=M{v\(x,t) 







We first determine the value of e such that, if ||t;o||ifi then classical solutions 
t;(-,i) to Navier-Stokes exist for all time. The argument holds for real t as well as 
in 



~Sy.^ {t:argte (-5,5)}, 



where < (5 < ^. Sectorial existence of analytic solution in t with exponential 
decay for large \t\ was useful in proving Theorem 12.21 We denote by At the class 
of functions analytic in t for t e Si for < |t| < T. 
We consider the space of functions 

where t = |i|e"^, \<j>\ < S, and the weighted norm 

\\v\\x= sup ||eJ^y-,t)||Hi+ sup j r||e3''y-,|i|e»^)||l,.d|4 . 
teSj,o<|t|<T \<p\<s [Jo " ) 

Note that 

/ \ / \ 



and / is the Fourier- Transform of /. 

The arguments below are an adaptation of classical arguments, see [3^. We 
introduce exponential weights in time, allowing for estimates independent of T, 
and extend the analysis to a complex sector. 

Lemma 12.2. For vq G H^, with zero average over T'^[0, 27r] we have 
where ci = ( 1 - 



Proof. First, take f = vo and t e [0,T]. Note that zero average imphes /(O) = 0; 

so we only need to consider > 1. 

(12.113) 

k^O k=£0 

Also, note that 

(12.114) £ \\e^'''e'''^f\\l.dt < J^l + \k\')\fkf (^j\-'''^^^^^'-i^'dt^ 

If |2< 2||.||2 
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li t ^ Sg, we integrate along the ray |i|e"^. It is clear all the steps go through when 
ly is replaced by i^cos0. A bound, uniform in Sg, is obtained by replacing ^ in 
(|12.114p by ^j^. The result follows. | 

Lemma 12.3. If ei^^F £ LhL^ uniformly in (j) £ (—6,6), then 



(12.115) 
with 

ci = 



f e'^(*-")^i^(x,T)dT 
Jo 

2V2 4V2 



<C2 sup ||e^*f^L2 

X U\<5 



yjv cos I i^cosSj 
Proof. We first show this for t E [0,T]. The function 

v{x,t) = f e''^'"''>'^F{x,T)dT 
Jo 

satisfies 

(12.116) vt-i'Av = F, v{x,0)=Q. 

Multiplying ()12.116p by v*, the conjugate of v, integrating over x G T^[0, 27r] and 
combining with the equation for v* we obtain 

(12.117) l|i^(.,i)||2,+2H|i?^^(-,i)iii. <^ll^^(-,<)lli. + ^lk(-,t)lli.. 

Similarly, taking the gradient in x of (|12.116p . taking the dot product with Vv* 
and combining with the equation satisfied by Vu*, we obtain 

^\\Dv{;t)\\h+MDM;mh= f {DF) ■ {Dv*)dx + f (DF*) ■ {Dv)dx. 
Integration by parts and Cauchy's inequality give 

(12.118) |pz;(.,<)||i.+2H|i?M-,<)lli^ < Jll^^(-,i)lli^ + ^l|At'(-,i)||i.. 
Combining p2.117p and p2.118p and using Poincare's inequality, we have 

(12.119) |ll<,t)lll^i+^HI<,i)llH.+ <^ll^(-,t)lli- 
Therefore, using (|12.119p and the fact that v{x, 0) = 0, 



vi;t)\\m<- I ||e^^n-,T)||i.dr. 



Hence, 



(12.120) sup ||ei'"t;(,i)|U. <^||e^'f^|L. ... 

tG[0,T] I*' 

Integration of (|12.119p . using u(x,0) = 0, gives 

rt 00 ft 



/ \\ei''^v{-,r)\\l.dT<% f |lei-n-,r)|li.dr. 
Jo ^ Jo 
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Therefore, for t e [0,T], we obtain 



(12.121) 



1/2 



4V2 



\\e-^'^^v{;T)\\l.Jr 
Now (|12.120p and p2.12ip together imply 

(12.122) sup |V(1 + Ifcp) e^-^'i / e-''l'=l'*(i^")/'(fc 



f\\ei''^F{;r)\\l.dr 
Jo " . 



1/2 



ts)ds 



2> 1/2 



d|t|5](l + |A:n 



< 



4 / ^(fc,is)e-'^''=' 
^2x/2 ^ 4^/2^^ '"^ 



2n 1/2 



1/2 



and replacing t e [0,T] by t e e*'^[0,T] £ is equivalent to replacing by lycosS. 
I 

Lemma 12.4. If F = -V [v ■ Wv], then for v eX, and t e e*'^[0,T] C Sg, 

sup llei'^'FL. <C3||H1^, 

I0I<5 " 



where C3 



/or i g S*?, and C4 is the Sobolev constant bounding \\ ■ ||^6 



by II ' ll-ffi T'See for instance P], pa^e 75j. 

Proof. First consider t G [0,r]. Holder's inequality implies 







1/2 




< 


r \e-'-^-\d\r\ 




[i 




JO 







|ei-|F(,r)|||i.d|r| 



1/2 



Hence, 



le^-^'Fll 



< 



(3t.)i/4 



lei-^'Fll 



If we replace t G [0,T] by t G 5*^ in this argument, the effect is simply that , K 



gets replaced by 



(3i.)i/4 



(3i/cos5)i/4 ■ 

For nonnegative u, w, repeated use of Holder's inequality gives 



vo^u^dx < 



T3 



< 



T3 



w^dx 



u^dx 



1/3 



1/3 



T3 



T3 



u^dx 



i^dx 



1/2 



2/3 



T3 



T3 



1/6 

dx \ < ||w|||6||m||l2||u||l5- 



Therefore, it follows that 



||e^*F(.,<)|U. < |lei''>(.,t)||V«(.,t)||U. < ||et''*HUsl|e^''*V^;||i//||et'^*V^;||i/^ 



and 



lei-^'FL. < |lef''*^;|U^..|le^*V^;||^/^^.||e^*Vt;||y.%,. 
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Using Sobolev inequalities, wc have 

< C4\\v{-,t)\\Hi, 
\\Dvi;t)\\Le<c4Dv{;t)\\Hl. 

Thus 

II ^i^tlPW ^ 3/2|| ivt ||3/2 I, 3 J 1/2 3/2|| ||2 

Therefore, 

*' (St^coso)^'* 

Since the right hand side is independent of taking the supremum of the left side 
over for |(/)| < 5, the Lemma follows. | 

Lemma 12.5. The operator M defined in 1112. 112\) satisfies the following estimate: 

\\^f[v]\\x < Cillwollffi +C2C3||w||x, 

||AA[«(i)] - AA[«(2)]||x < C2C3 (ik^^^lU + Wv^'^Wx) \\v^'^ - v'^'^'Wx- 
Proof. Note that 

AA[i;]=e''*^wo+ / e'^(*-")^F(-,T)dT, 
|/o 

where F = -V [v ■ Vv]. By Lemmas [1221 [1231 and [123 it follows that 

IIA^H \\x < Cill^olUi +C2C3||w|||. 

For the second part, we note that 

Using Lemmas 112.21 [lX51 and [T2.4I again, we obtain the desired estimate. | 
Lemma 12.6. // 

1 3l/4l.7/4[cos5]7/4 



4ciC2C3 8V2c4/^(\/j/cos5 + V2)(2 + \/z/cos5)' 



v{x,t) exists in X for any T. v{-,t) is analytic int d Sg and decays exponentially 
in that sector as \t\ — > oo, with 

Further, this solution is smooth in x. If 

3l/4i.7/4 

Ikollffi < eo = 



8V2cf'(V^+V2)(2 + V^)' 
then v(x, t) is a classical solution for all t G K+. 

Proof. If ||fo||i/i < e, Lemma 112.51 implies that the operator A/" (defined in Lemma 
112. 5p is contractive and hence a solution to Navier-Stokes equation exists in X. 
Since the estimates are uniform in <, it follows that this solution exists for all 
t G Sg. Known results (or Theorem 12.11 above) imply that if the initial data is in 
H^, then the solution becomes smooth (in fact, analytic for periodic data, [TB]) 
instantly, and thus it is a classical solution when t > 0. Analyticity and decay 
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in t follows from the definition of X, the arbitrariness in the choice of T and the 
observation that A/" in Lemma 112.51 is contractive in a ball of radius 2ci||uo||ffi- 
Further, by taking the lim^^Q+ e = eg, we obtain the less restrictive condition on 
||wo|l_f/i that ensures existence of classical solution only for t S K"*". | 

Lemma 12.7. // HwoHh^ < £2 for sufficiently small 62, 
for any t €L Sg. 

Proof. This is similar to the proof of Lemma 112.61 with X replaced by 
X = {Aif'}n{Lf,|i73| (^At®H^{T^[0,2n])^ n(^L^ [e'^iO,T)]®H^{T^[0,2TT] 

for 101 < ^. I 

Theorem 12.1. A weak solution to U.l]) becomes classical when t > T^., where 

_ 256^ct(VI7+\/2)2(2 + yi7)2 



31/2^,9/2 

This solution is analytic in t for (t — Tc^a) G Sg, where 

^ _ 256Ecl(Vi^cosS + a/2)2(2 + VTcosS)^ 

~ 3l/2z.[j.c0S^]7/2 ■ 

Further, for any constant C , there exists T2 so that for {t — T2) G Sg, 

\\v{-Mi^ < Cexp 
Proof. Leray's energy estimate implies 



\Re{t-T2} 



where E = ^\\vo\\^^2- From a standard pigeon-hole argument, it follows that there 
exists Ti e (0, T] so that 

||V<,Ti)|U. <yj. 
Therefore, Poincare's inequality implies 

This means there exists some Ti G [0,Tc], where 



31/2^,9/2 
for which 

31/4^7/4 

""^•''^^^""^%V2cr(VI^+V2)(2 + yi7)- 
Replacing thy t — Ti in Lemma fl 2. 61 we see that the solution is classical and smooth 
for t — Ti E K.+ , therefore necessarily for t > Tc. 
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Further, from these arguments, it is clear that there exists a Ti^a G [0,Tc,a] so 
that 

n . M, 3l/4[i. COS ^17/4 

8V2cl'^ViycosS + V2){2+ ViycosS) 
Replacing t by t — Ti ^ in Lemma (|12.6I) . we see that the classical solution is analytic 
in t — Ti c G Sg (which includes the region t — Tc^a G Sg). 
Further, since for t > Ti we have 

\\v{-,t)\\],,dt< sup ll^-lll < (2cieo)2, 

Ti T>Ti 

it follows from a pigeon-hole argument that given 62, there exists a T2 > Ti such 
that 

\\v{-,T2)\\h^ < £2- 

From Lemma 112.71 it follows that v exists for t — T2 E Sg and 

||t;(.,t)||^. <2cie2e~^«'=(*-^^). 
The last part of the theorem follows from (recall ^(O) ~ 0) 

\\H-,t)y <c4\k\M-,t)y <c5\\v{-,mm- 

I 

Remark 12.8. The decay rate e~i'^* for ||w(-,i)|j/i is not sharp. A more refined 
argument can be given, to estimate away the nonlinear terms and obtain a e~^* 
decay. 
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Figure 2. For zero forcing and f = 0.1: (a). The original (unac- 
celerated) equation, vs. p. (b). Accelerated equa- 

tion with n = 2, vs. q. 



NAVIER-STOKES EQUATION 
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where s = q^^^ and A_ is the backward difference operator in s. 
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